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FINE AND HYPERFINE STRUCTURE 
OF THE MUONIC 3 He ION 

E.N. Elekina, A.P. Martynenko* 

Samara State University, Pavlov Street 1, Samara 443011, Russia 

On the basis of quasipotential approach to the bound state problem in QED we 
calculate the vacuum polarization, relativistic, recoil, structure corrections of or- 
ders a 5 and a 6 to the fine structure interval AE^ S = i?(2P 3 / 2 ) — £ , (2P 1/ / 2 ) and to 
C^-r the hyperfine structure of the energy levels 2P\/2 and 2P 3/ / 2 in muonic \He ion. 

The resulting values AE* 8 = 144803.15 (ieV, AE h f s (2P l/2 ) = -58712.90 fxeV, 
(y-^ . AE^ s (2P 3 / 2 ) = —24290.69 fieV provide reliable guidelines in performing a com- 

parison with the relevant experimental data. 

PACS numbers: 31.30.Jv, 12.20.Ds, 32.10.Fn 



in: 



in 

(N 



o 



Keywords: fine and hyperfine structure, muonic helium ion 



INTRODUCTION 



Simple atoms play important role in the check of quantum electrodynamics (QED), the 
bound state theory and precise determination of fundamental physical constants (the fine 
structure constant, the lepton and proton masses, the Rydberg constant, the proton charge 
radius, etc) [H, 0, y| . Light muonic atoms (muonic hydrogen (fip), muonic deuterium, ions 
of muonic helium etc.) are distinguished among simple atoms by the strong influence of the 
vacuum polarization (VP) effects, recoil effects, nuclear structure and polarizability effects 
on the structure of the energy levels. The comparison of the theoretical value of the fine and 
hyperfine splittings in muonic helium ions with the future experimental data will lead to a 
more precise value of the helion charge radius and the check of quantum electrodynamics 
with the accuracy 10~ 7 . The energy levels of muonic helium ions were theoretically studied 
many years ago in jl, 0, @, 0, 0] both on the basis of the relativistic Dirac equation and 
nonrelativistic approach, accounting different corrections by the perturbation theory (PT). 
In these papers the basic contributions to the energies for the (2P — 2S) transitions in muonic 
helium (n\He) + were evaluated with the accuracy 0.1 meV. 

In this work we continue the investigation [9| of the energy spectrum of {^\He) + in the P- 
wave part. The aim of the present study is to calculate such contributions of orders a 5 and a 6 
both in the fine and hyperfine structure of the energy states 2Pu2, 2P3/2, which are connected 
with the electron vacuum polarization, the recoil and structure effects, the muon anomalous 
magnetic moment and the relativistic corrections. The role of all these effects is crucial in 
order to obtain high theoretical accuracy. Our purpose also consists in the refinement of the 
earlier performed calculations in 0, M, @, S] and in the derivation of the reliable numerical 
estimate for the structure of P-wave levels in the ion {^\He) + , which can be used for the 
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comparison with experimental data. Modern numerical values of fundamental physical con- 
stants are taken from Ref.j^|: the electron mass m e = 0.510998910(13) • 10~ 3 GeV, the muon 
mass m^ = 0.1056583668(38) GeV, the fine structure constant a' 1 = 137.035999679(94), 
the proton mass m p = 0.938272013(23) GeV, the helion mass 2.808391383 (70) GeV, the 
helion magnetic moment \Xh = —2.127497723(25), the muon anomalous magnetic moment 
= 1.16592069(60) • 10~ 3 . 



II. FINE STRUCTURE OF P - WAVE ENERGY LEVELS 



Our approach to the investigation of the energy spectrum of muonic helium ion {ix\He) + is 
based on the use of quasipotential method in quantum electrodynamics [l(J 11, 121 ]. where the 
two-particle bound state is described by the Schrodinger equation. The basic contribution 
to the muon and proton interaction operator is determined by the Breit Hamiltonian 13, 
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where m±, mi are the muon and proton masses, fi 
AV fs is the muon spin-orbit interaction: 
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AV h f s is the helion spin-orbit interaction and the interaction of the muon and helion spins. 
The leading order {Za) 4 contribution to the fine structure is determined by the operator 
AV^ S . As it follows from Eq.(2), the potential AV^ S includes also the recoil effects (the 



Barker-Glover correction [16]) and the muon anomalous magnetic moment correction. 
The fine structure interval (2P 3 / 2 — 2Pi/ 2 ) for the ion (/^ife)" 1- can be written in the form: 
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This expression includes a relativistic correction of order {Za) G , which can be calculated 
with the aid of the Dirac equation [j], [13], the correction of order a{Za) 6 enhanced by 
the factor \n{Za) [H, 19, 2^, a number of terms of fifth and sixth order in a which are 
determined by the effects of the vacuum polarization. The relativistic recoil effects of order 
rju{Zaf jm 2 in the energy spectra of hydrogenic atoms were investigated in Refs. [3, E3, 21 



1221 . 1231 ] . In the fine splitting (3) they were calculated in [171 . 123] . Additional corrections of 
the same order were obtained in [24]]. They do not depend on the muon total momentum 
j and give the contribution only to the Lamb shift. The contributions to the coefficients 
Ayp and Byp arise in the first and second orders of perturbation theory. Numerical values 
of the terms in the expression (3), which are presented in the analytical form, are quoted 
in Table I for definiteness with the accuracy 0.01 [ieV. The fine structure interval (3) in 
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FIG. 1: One-loop vacuum polarization contributions to the fine and hyperfine structure. The 
dashed line corresponds to the Coulomb interaction. The wave line corresponds to the fine or 
hyperfine interaction. G is the reduced Coulomb Green's function. 



the energy spectrum of electronic hydrogen is considered for a long time as a basic test of 
quantum electrodynamics [13, 25. 2a . 

The leading order vacuum polarization potential which gives the contribution to the 
coefficient Ayp, is presented by the Feynman diagrams in Fig.l. The one- loop vacuum 
polarization effects lead to the modification both the Coulomb interaction and the spin- 
orbit interaction in expressions (1), (2) [l3|, 14j : 
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where the spectral function p(s) = \/s 2 — l(2s 2 + l)/s 4 , m e is the electron mass. Averaging 
the potential (2) over the wave functions of the 2P - state 
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we obtain the following contribution to the interval (3) (see Fig. 1(a)): 
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Although the integral in Eq.(7) can be calculated analytically, we present here for simplicity 
only its numerical value. 

Higher order corrections a 2 (Za) 4 entering in the are taken into account in this ex- 
pression as well as the recoil effects. The same order contribution a(Za) 4 can be obtained 
in the second order perturbation theory (see Fig. 1(b)). In this case the energy spectrum is 



determined by the reduced Coulomb Green's function [14|, |27|, [28 
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where z < = min(z, z'), z> = max(z, z'), C = 0.577216... is the Euler constant, z = Wr. 
Using Eqs. (8) and (9) we transform the correction of order a(Za) A to the fine structure in 
the second order perturbation theory as follows: 
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Let us consider the two-loop vacuum polarization contributions in the one-photon inter- 
action shown in Fig. 2. They give the corrections to the fine splitting of P- levels of order 
a 2 (Za) 4 . 




FIG. 2: Effects of two-loop electron vacuum polarization in the one-photon interaction. 



In order to obtain the particle-interaction operator for the amplitude, corresponding to 
the diagram in Fig. 2 (a), it is necessary to make the substitution 
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two times in the photon propagator. In the coordinate representation, the interaction oper- 
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Averaging (12) over wave functions (6), we obtain the correction to the interval (3): 
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The two-loop vacuum polarization operator is needed to find the interaction operator shown 



in Fig.2(b,c). The modification of the photon propagator in this case has the form [31 
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The two-loop vacuum polarization potential and the correction to the fine structure (2P3/2 
2P1/2) are the following: 
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FIG. 3: Effects of two-loop electron vacuum polarization in the second order perturbation theory. 
The dashed line corresponds to the Coulomb interaction. The wave line corresponds to the fine or 
hyperfine interaction. G is the reduced Coulomb Green's function. 

Two-loop vacuum polarization contributions in the second order perturbation theory 
shown in Fig. 3, have the same order a 2 (Za) 4 . In order to calculate them, it is necessary 
to employ relations (2b (J^ (5), (8), and the modified Coulomb potential by the two- loop 
vacuum polarization [lOj, 
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The amplitude (a) in Fig.3 gives the following correction of order a 2 (Za) 4 to the fine split- 
ting: 
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Two other contributions from the amplitudes (b), (c) in Fig.3 have the similar integral 
structure. Their numerical values are included in Table I. The summary result for the fine 
splitting AE^ S in (fi 2 He) + is presented here also. It takes into consideration the numerous 
earlier performed calculations discussed in the review article [l| and new corrections obtained 
in this work. 



III. HYPERFINE STRUCTURE OF THE ENERGY LEVELS 2P 1/2 AND 2P 3/2 



The leading order contribution to the hyperfine splitting of the energy levels 2Py 2 and 
2P 3 / 2 in muonic helium ion (u^He) + of order a 4 is determined by the potential (the hyperfine 
part of the Breit potential) [131 ] : 
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where n = r/r. The operator (21) does not commute with the operator of the muon total 



angular momentum J 



\cr v This leads to the mixing of the 2P\j 2 and 2P 3 / 2 energy 



levels and, hence, to a more complicated hyperfine structure of P-wave levels. 

In order to calculate the diagonal matrix elements (2P1/2I AVg |2Pi/ 2 ) and 
(2P3/2I AV^^ s |2P 3 /2) we can use the following replacements for the operators (siS 2 ) and 



(Ls 2 ), which involve the spin of the nucleus [3 
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where (si«J), (LJ) are eigenvalues of corresponding operators between the states with equal 
orbital momentum /. In addition, the averaging over angles in the second term in the 



right-hand side of (21) can be carried out by means of the relation [32 
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The diagonal matrix elements (2P 1/2 | AV hfs \2P 1/2 ) and (2P 3/2 | AV hfs \2P 3/2 ) lead to the fol- 
lowing hyperfine structure: 

AP^(2P 1/2 ) = P(2 3 P 1/2 ) - E(2 1 P 1/2 ) = (24) 
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TABLE I: Fine structure of P-wave energy levels in muonic \Ue ion. 
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where Ep = Z 3 a 4 ii 3 iih/3mim p is the Fermi energy for the n = 2 level. The calculation of 
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the relativistic corrections A^J, A^J within this approach includes the study of two-photon, 
three-photon exchange diagrams and the second order perturbation theory contributions 
with the Breit Hamiltonian determined by Eqs. (1), (2) and (21). More simple method 
for their calculation is based on the relativistic Dirac equation [20 
hyperfine interaction potential has the form: 



331 ] . In that case, the 
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where the nuclear magnetic moment fi = Qn^n s 2 (a*jv = e/2m p ). The radial integrals 
Rk — Io° dkfkdr, which are characteristic of this case, are determined by the Dirac wave 
functions fk,gk of the states 2Pi/ 2 and 2P 3 / 2 . Taking into account their explicit form 32 



we obtain the following values of the relativistic corrections to the hyperfine structure of the 
P - wave levels: 
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These values of the coefficients coincide with the result obtained in Ref. |23j by analytically 
calculating the contribution of order m\(Za) Q /ni2 to the hyperfine structure of the P - wave 
levels of the hydrogen atom at n = 2. 

The fifth order contribution over a due to the electron vacuum polarization (see diagrams 
(a), (b) in Fig.l) appears in the hyperfine splitting in just the same way as in the fine 
structure of the spectrum. The modification of the hyperfine part of the Breit potential 
induced by the vacuum polarization is described by the following expression (the substitution 
(11) is used) [l7 
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The subsequent transformations of the diagonal matrix elements of the operator (29), con- 
nected with the averaging over angles, can be performed with the use of the formula (23) 
for the second term in the square brackets of Eq.(29). Similar averaging for the first term 
in the square brackets is: 

(5ij - riiTij) = -(25ij + LiLj + LjLi). (30) 

Then the contributions of the vacuum polarization in the first and second orders of pertur- 
bation theory can be written as follows: 
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Two-loop vacuum polarization corrections to the hyperfine part of the potential for the 
/ 7^ can be obtained with the aid of relations (11) and (14). The results are 



Av hfs 



Za/jf, 



vp-vp\ 



2m 1 m 2 r 



a\ 2 r°° r°° 1 

-) / Mdtl PW,—^ 



(35) 



1 + 



mi Zmim v 



m 2 



2m 2 2 p h 



(LtTr 



£ 2 (1 + 2m e £r)e~ 2me€r - r] 2 (l + 2m e r/r)e~ 



2m e rir 



1 + a. 



CT \ CT 1 



3(o- 1 n)(cr 2 n)) (£ 2 (1 + 2m e £r)e~ 2m ^ r - r? 2 (l + 2m e T 1 r)e 



■2m e r]r 



+ 



+4m 2 r 2 {<j x <j 2 - S^nX^n)) (^V 2 ™^ - v 4 e~ 2m ^ 



W 2 h lL vyP {r) 



Za^h 2 fa 



2m,\m 2 r z 3 \7r 



1 f(v)dv -^psz= 
■ e x 



x • 



1 + 



mi Zmiirir 



m 2 2m\p h 



1 + 



2m P r 



o 1 — 
(Lcr ~ 



(36) 



Am 2 r 2 . 
1^ 



- 1 n)( < r 2 n))+ 1 + 



2m P r 



( (T i (T 2 



1 + a, 



3(< T i n )( cr 2 n )) 



Vl^v 2 , 

Omitting further details of the calculations that can be performed with the aid of a procedure 
similar to that which was used in the derivation (17) and (20), we represent in Table II the 
numerical values of the contributions to the energy spectrum that are determined by the 
potentials (35) and (36). Yet another part of two-loop corrections to the hyperfine structure 
in the second order PT is shown in Fig. 3. We also included in Table II numerical results of 
the contributions from these amplitudes to hyperfine structure of P-wave states. 

Nuclear structure effects play significant role in the precise calculation of the hyperfine 
structure in muonic atoms [9j, 10 . In order to calculate it in the case of P-wave levels we 



can expand the helion magnetic form factor over relative momenta and obtain the hyperfine 
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TABLE II: Hyperfine structure of P-wave energy levels in muonic helium ion (n 2 He) + . 
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-58356.61 


-24088.50 
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Muon AMM contribution 
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VP Contribution in the second 


-0.45 


-0.19 


(11),(19),(21) 


order PT of order a 6 








(AV 2 c _ loopyP -G-AV h f°) 








Nuclear structure correction 


-0.33 


0.66 


(38) 


Summary contribution 


-58539.90 


-24117.69 





part of the interaction operator in the momentum representation, which is proportional to 
magnetic mean-square radius < r 2 M > of the helion: 



AKl /s (k) 



wZa(l + a M ) < r M > 
6mi?Ti2 



> lC r 2 )k 2 - (<r 1 k)(o- 2 k) 



(37) 



Averaging this operator over the wave functions ip2P we obtain the following contribution 
to the hyperfine structure: 



AE 



hfs 
sir 



/i 5 (Za) 6 (l + a„) < rl > [F(F + 1) - j(j + 1) - + 1) - f ] 



144mim2 



(38) 



Corresponding numerical values for the hyperfine splittings of 2Pi/ 2 and 2P3/2 states are 
included in Table II. We assume that the values of the charge and magnetic mean-square 
radii coincide and take the value of 2 He charge radius te = 1.9642(11) fm 34 . 
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2P 



3/2 



2P 



1/2 



2S 



1/2 



A 




AE is 



AE 



Ls 



2 3 P 



3/2 



AP hfs (2P 3/2 ) 
2 5 P 



3/2 




2 1 R 



1/2 



AP hfs (2P 1/2 ) 



2 3 P 



1/2 




2^ 



1/2 



AP hfs (2S) 



2 3 S 



1/2 



FIG. 4: The structure of 5-wave and P-wave energy levels in muonic helium ion (/x 2 Pe) + for the 
n = 2. 



The off-diagonal matrix element has an important role to attain the high accuracy of the 
calculation of the P-wave levels in muonic helium ion. We present its general structure as 
follows: 



2m\ Zmirrir, 



+ -fC rel (Za) 2 + C VP a 



7 = <2 3 P 1/2 |AV^|2 3 P 3/2 > = E F {-^ \l - a, ^ 

(39) 

where for the sake of simplicity, we have restricted ourselves to terms of fifth order in a 
in considering vacuum polarization effects, terms of fifth and higher orders in the muon 
anomalous magnetic moment, relativistic effects of order {Zotf and recoil effects. The first 
three terms in the right-hand side of Eq. (37) result from employing the potential (21). In 
the Dirac's theory, relativistic corrections are determined by off-diagonal radial integrals: 



(40) 
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By using the explicit expressions for the wave functions fi/ 2i 3/ 2 (r) and #1/2,3/2 ( r ) [32J in 
order to calculate this integral, we obtain the coefficient C re \ = 9/16. In order to obtain the 
vacuum polarization correction in 7 we use the potential (29). Then we have to calculate 
the matrix elements of the following operators: 



T-i 



(L<r 2 ), T 2 



[ (T l (T 2 



3(o- 1 n)(<r 2 n)] , T 3 



[ (T l (T 2 



V Cr i n )( Cr 2 n )] 



(41) 



Upon averaging over angles with the aid of expressions (23) and (30), they can be written 
in terms of the 6j - symbols as 



(T 3 ) = -(T 2 ) = -(T 1 )== 6jj'< 



IF 1 

.22-7 



1 F 1 

2 2-7. 



2V2 



(42) 



where the value of the total momentum F = 1 (F = s 2 + J), I = 1, j = v/2.7 + 1, f = 
\ / 2j' + 1, and the numerical values of the 6j - symbols are borrowed from [35(. As a result 
the vacuum polarization contributions to the off-diagonal matrix element (39) in the first 
and second orders of perturbation theory have the form: 



7i 



(2 3 Pi 



1/2 



AVy p s \2 3 P 3/2 ) = E F 



P(s)ds 



xdxe 



-x( 1- 



x 



1 + ^ 

m 2 



Zminip 
2m\ii h 



2m P s 
1 + -r^x 



1 + a u ( 2m e s Amis 2 .-> 

r — 1 + -r^X t^X 1 



72 = (2 3 P 



1/2 



AV^ p -G-AVb Is ^ 3 



x 



p(s)ds 



W 

2 a P 3 /2> = E F (^ 

(-1 . 2m e s \ 




W W 2 



V2\ 



a 



2592 J 7T 
dx' 



2m\ Zmirrip 



m 2 



rn£nh 



(43) 

9.49 fjteV, 
x (44) 



x 



g(x,x 



5.33 neV. 



The total numerical value of the matrix element (39) is 7 = 5497.28 fieV. It leads to the shift 
of the hyperfine splittings of the energy levels 2 3 P 3 / 2 and 2 3 P!/ 2 by the value 5 = 173.0 fieV. 



IV. SUMMARY AND CONCLUSION 

In the present study we have calculated QED effects in the fine and hyperfine structure 
of the 2P1/2 and 2P 3 / 2 energy levels in muonic helium ion {p\He) + . We have considered the 
electron vacuum polarization contributions of orders a 5 , a 6 , recoil corrections, relativistic 
effects of order a 6 and nuclear structure corrections. The numerical values of the contribu- 
tions are presented in Tables I and II. In these Tables we give the references to other papers 
also devoted to the investigation of the structure of P-wave levels in the hydrogenic atoms. 

Let us summarize the basic particularities of the calculation performed above. 

1. Special attention in our investigation has been concentrated on the vacuum polar- 
ization effects. For this purpose we obtain the terms of the interaction operator in muonic 
helium ion which contain the one-loop and two-loop vacuum polarization corrections. 

2. In each order in a we have taken into account recoil effects in the terms proportional 
to mi/m 2 . The experimental values of the muon and helion magnetic moments are used 0]. 

3. The calculation of the relativistic corrections to the diagonal and nondiagonal matrix 
elements both for the fine and hyperfine structure intervals is performed on the basis of the 
Dirac equation. 
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TABLE III: (25 — 2P) transition energies for muonic helium ion ((j, 2 He) 



Transition 


Energy (meV) 


M 


1 S\/2 — 3 Pl/2 


1167.33 


1167.3 


1 S\/2 — 3 P3/2 


1342.19 


1342.3 


3 5'l/2 — 1 Pl/2 


1392.66 


1392.4 


-J 1/2 ~~ -^1/2 


1333.95 


1334.1 


3 5"i/2 — 3 P3/2 


1508.80 


1509.1 


3 Sl/2 — 5 -P3/2 


1484.51 


1484.6 



Total numerical values for the fine structure interval AE? S (3) and hyperfme splitting 
intervals of 2Pi/2 and 2P3/2 states are presented in Tables I, II. Taking also into account our 
calculation of the mixing of the 2 3 P 3 / 2 — and 2 3 Pi/ 2 — wave energy levels (the correction 5), 
we find that these values of hyperfme splittings change by 5 = 173.0 fieV: AE h ^ s (2Pi/ 2 ) = 
AE hfs (2P 1/2 ) - 5 = -58712.90 ^eV, AE hfs (2P 3/2 ) = AE hfs (2P 3/2 ) - 5 = -24290.69 //eV '. 
The theoretical error of the obtained results is determined by the contributions of higher 
order and amounts up to 10 -6 . Previously, the (2S — 2P) transition energies for muonic 
helium ion were studied in 0, 0, 0]. Considering the obtained in 0, 0] numerical results 
for different transitions (2s — 2p) we find that the fine splitting and the hyperfme split- 
ting intervals for the states 2P\/ 2 and 2P 3 / 2 in this paper are equal: AE? S = 145.0 meV, 
AE hfs (2P 1/2 ) = 58.3 meV, AE hfs (2P 3/2 ) = 24.5 meV. So, the results of our work agree 
and refine the previous calculations performed in p, 0] via taking into account higher order 
effects. They can be considered as a reliable estimate for the fine and hyperfme structure 
intervals for the P- levels in muonic helium ion (n\He) + . The disposition of the P-wave 
energy levels is shown in Fig.4. Taking into account the value of the Lamb shift (2^ — 2P) 
from our result for the hyperfme splitting of 2S state from [9| and the numerical results 
obtained in this work, we find new (2S* — 2P) transition energies, which are presented in 
Table III. 
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